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Abstract. Harmonic weak Maass forms of half-integral weight are the subject of many 
recent works. They are closely related to Ramanujan's mock theta functions, their theta 
lifts give rise to Arakelov Green functions, and their coefficients are often related to central 
values and derivatives of Hecke L- functions. We present an algorithm to compute harmonic 
weak Maass forms numerically, based on the automorphy method due to Hejhal and Stark. 
As explicit examples we consider harmonic weak Maass forms of weight 1/2 associated to 
the elliptic curves llal, 37al, 37bl. We made extensive numerical computations and the 
data we obtained is presented in the final section of the paper. We expect that experiments 
based on our data will lead to a better understanding of the arithmetic properties of the 
Fourier coefficients. 



1. Introduction 

Half-integral weight modular forms play important roles in arithmetic geometry and 
number theory. Their coefficients serve as generating functions for various interesting 
number theoretic functions, such as representation numbers of quadratic forms in an odd 
number of variables or class numbers of imaginary quadratic fields. Moreover, employing 
the Shimura correspondence [Shj . Waldspurger |Wa] . and Kohnen and Zagier \KZ\ IK] 
showed that the coefficients of half-integral weight cusp forms essentially are square-roots 
of central values of quadratic twists of modular L-functions. In analogy with these works, 
Katok and Sarnak [KSJ used a Shimura correspondence to relate coefficients of weight 1/2 
Maass forms to sums of values and sums of line integrals of Maass cusp forms. 

In more recent work Zagier discovered that the generating function for the traces of 
singular moduli (the CM values of the classical j-function) is a weakly holomorphic modular 
form of weight 3/2 |Zalj . This result, which was generalized in various directions (see e.g. 
|B02j . |BF2j . |DJj . |Kij ) . demonstrates that also the coefficients of automorphic forms with 
singularities at the cusps carry interesting arithmetic information. 

In a similar spirit, Ono and the first author proved that the coefficients of harmonic weak 
Maass forms of weight 1/2 are related to both the values and central derivatives of quadratic 
twists of weight 2 modular L-functions [BruO] . Harmonic weak Maass forms are also closely 
related to mock modular forms and to Ramanujan's mock theta functions, which have been 
the subject of various recent works (see e.g. |BUll \B03\ lOnl IZa2l \Zwl\ IZw2j ). In view of 
these connections, it is desirable to develop tools for the computation of such automorphic 
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forms. In the present paper we propose an approach to this problem which yields an 
efficient algorithm. Moreover, we compute some harmonic weak Maass forms which are 
related to rational elliptic curves as in |BruO] . 

The non-holomorphic nature of harmonic weak Maass forms prevents the use of the 
well developed algorithms existing for (weakly) holomorphic modular forms, such as e.g. 
modular symbols. The use of Poincare series does not work well either in small weights 
due to the poor convergence of the infinite series which appear in the explicit formulas 
for the coefficients. Instead we adapt the 'automorphy method', originally developed by 
Hejhal for the computation of Maass cusp forms on Hecke triangle groups (see e.g. |He] ). 
to the setting of harmonic weak Maass forms. 

We now describe the content of this paper in more detail. Let k G and let be a 
positive integer (with 4 | if G |Z \ Z). A harmonic weak Maass form of weight k on 
ro(A^) is a smooth function on H, the upper half of the complex plane, which satisfies: 

(i) /U7 = /forall7ero(Ar); 

(ii) Afc/ = 0, where is the weight k hyperbolic Laplacian on H (see (12. 3p ): 

(iii) There is a polynomial Pj = J2n<o^^i^)^'^ ^ such that /(r) — Pjir) = 
0{e~'^^) as f — )■ oo for some e > 0. Analogous conditions are required at all cusps. 

Throughout, for r G H, we let r = m + iv, where M,f G M, and we let q := e^'^*'^. The 
polynomial Pj is called the principal part of / at oo. 

Such a harmonic weak Maass form / has a Fourier expansion at infinity of the form 

(1.1) /(r)= c+(n)g" + 5^c-(n)r(l-fc,47r|n|t;)g", 

n3>— oo n<0 

where T{a,x) denotes the incomplete Gamma function. The series ^„^_oo 

called the holomorphic part of /, and its complement is called the non-holomorphic part. 

Naturally, / has similar expansions at the other cusps. There is an antilinear differential 

operator, taking / to the cusp form ^k{f) '■= 22u*^|| of weight 2 — k, see (12. 5p . The kernel of 
consists of the space of weakly holomorphic modular forms, those meromorphic modular 
forms whose poles (if any) are supported at cusps. 

Every weight 2 — k cusp form is the image under C,k of a weight k harmonic weak Maass 
form. Ramanujan's mock theta functions correspond to those forms whose images under 
.^1/2 are weight 3/2 unary theta functions. Here we mainly consider those weight 1/2 
harmonic weak Maass forms whose images under C,i/2 are orthogonal to the unary theta 
series. According to |BruO] . their coefficients are related to both the values and central 
derivatives of quadratic twists of weight 2 modular L-functions. 

We now briefly describe this result in the special case that the level is a prime p. Let G G 
S2{Tq{p)) be a normalized Hecke eigenform whose Hecke L-function L{G, s) satisfies an odd 
functional equation. That is, the completed L-function A(G,s) = p*/^(27r)~*r(s)L(G, s) 
satisfies A(G, 2 — s) = Eg-^^G, s) with root number eq = —1- Therefore, the central critical 
value L{G,1) vanishes. By Kohnen's theory of plus-spaces [Kj, there is a half-integral 
weight newform g G S^^^{To{4:p)), unique up to a multiplicative constant, which lifts to 
G under the Shimura correspondence. We choose g so that its coefficients are in Fq, the 



COMPUTATION OF HARMONIC WEAK MAASS FORMS 



3 



totally real number field generated by the Hecke eigenvalues of G. There exists a weight 
1/2 harmonic weak Maass form / on ro(4p) in the plus space whose principal part Pf has 
coefficients in Fg-, and such that 

6/2(/) = ikir^^, 

where \\g\\ denotes the usual Petersson norm. For a fundamental discriminant A let xa be 
the Kronecker character for Q(a/A), and let L{G, xa, s) be the quadratic twist of L{G, s) 
by Xa- One can show that the root number of L{G, xa, s) is equal to sign(A) ■ Xa(p) ^g- 

Theorem 1.1 (See [BruOj ). Assume that G, g, and f are as above, and let c^{n) denote 
the Fourier coefficients as in fll.ip . 

(1) //A < is a fundamental discriminant for which = 1, then 

L{G,XAA)=8n'\\Gr\\gr^-c-{Ar. 

(2) // A > is a fundamental discriminant for which = 1, then L'(G,xa, 1) = 
if and only z/c^(A) is algebraic. 

Note that the harmonic weak Maass form / is uniquely determined up to the addition of 
a weight 1/2 weakly holomorphic modular form with coefficients in Fq- Furthermore, the 
absolute values of the nonvanishing coefficients c"'"(A) are typically asymptotic to subex- 
ponential functions in n. For these reasons, the connection between L'{G,xa, 1) and the 
coefficients c^(A) in Theorem 11.1( 2) cannot be modified in a simple way to obtain a for- 
mula as in the first part of the Theorem. In fact, the proof of Theorem 11.1( 2) is rather 
indirect. It relies on the Gross-Zagier formula and on transcendence results of Waldschmidt 
and Scholl on periods of differentials on algebraic curves. 

The above result is one of the main motivations for the present paper. Our goal is to carry 
out numerical computations for the involved harmonic weak Maass forms. In that way we 
hope to find more direct connections of the coefficients c^(A) to periods or L-functions. 
When L'{G, Xa, 1) vanishes, meaning that c^(A) is algebraic (actually contained in Fg), it 
would be interesting to see if c"^(A) carries any arithmetic information related to G. In a 
forthcoming paper |Br2] . the coefficients c'^{n) will be linked to periods of certain algebraic 
differentials of the third kind on modular curves. It leads to a conjecture on differentials of 
the third kind on elliptic curves, which is based on the numerical data presented in Section 
IHof the present paper. 

Our computations make use of an adaption of the so-called automorphy method. The 
key point of this method is to view an automorphic form on a non-co-compact (but co- 
finite) Fuchsian group F as a function on the upper half-plane with certain transformation 
properties under the group F as well as convergent Fourier series expansions at all cusps. 
This classical point of view, in terms of functions on the upper half-plane, stands in con- 
trast to the more algebraic point of view, in terms of Hecke modules, usually taken when 
computing holomorphic modular forms. 
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By computing an automorphic form in this setting we mean that to any given (small) 
e > we compute a sufficient number of Fourier coefficients, each to high enough precision, 
so that we are able to evaluate the function at any point in the upper half-plane with an 
error at most e. 

To calculate these Fourier coefficients we truncate the Fourier series representing and 
view the resulting trigonometric sum as a finite Fourier series. Using the Fourier inversion 
theorem together with the automorphic properties of (which will additionally intertwine 
the Fourier series at various cusps) we are able to obtain a set of linear equations satisfied 
approximately by the coefficients. Cf. e.g. |Het IStlt IAv2] . The (surprising) effectiveness of 
this algorithm is closely related to the equidistribution properties of closed horocycles (cf. 
e.g. |Helt IB]). We describe the main algorithm in detail in Sectional The implementation 
of the software package is briefly described in Section 13.31 

In Section m we describe our computational result in three cases of particular interest. 
We consider the elliptic curves llal, 37al, and 37bl and their corresponding weight 2 
newforms. For instance, the elliptic curve 37al, is the curve of smallest conductor with 
rank 1. It corresponds to the unique weight two normalized newform G on ro(37) whose L- 
function has an odd functional equation. We verified the statement of Theorem II. II for all 
fundamental discriminants A which are squares modulo 148 in the range < A < 15000. 
For eight of these fundamental discriminants the quantity L'{G, xa, 1) vanishes. In all these 
cases we found a stronger statement then that of the Theorem 1 1.1 1 to be true, namely, that 
the associated coefficient c"'"(A) was an integer. For the corresponding data see Tables HI 
and [5l We conclude Section |4| by describing some analogous experiments for newforms G 
of weight 4, where g is of weight 5/2 and / of weight —1/2. 

The present paper is organized as follows. In Section [2] we recall some facts on (half 
integral weight) harmonic weak Maass forms. When working with arbitrary (not necessarily 
prime) level, it is convenient to use vector valued modular forms. In Section I2l3] we therefore 
recall from [BruOj the vector valued version of Theorem 11.11 In Section [3] we describe the 
automorphy method in the context of harmonic weak Maass forms. In Section HI we collect 
our computational results. In particular, we present results for the elliptic curves llal, 
37al, and 37bl; cf., e.g. Tables [H S] and [71 More extensive tables can be obtained from the 
authors on request. 

2. Preliminaries 

In order to be able to work with newforms of arbitrary level, it is convenient to work 
with vector valued modular forms of half integral weight for the metaplectic extension of 
SL2(Z). We describe the necessary background in this section. 

2.1. A Weil representation. Let H = {r G C; Q'(t) > 0} be the complex upper half 
plane. We write Mp2(M) for the metaplectic two- fold cover of SL2(M), realized as the group 
of pairs (M, 0(r)), where M = (° ^) G SL2(M) and : H — C is a holomorphic function 
with 0(t)^ = ct + d. The multiplication is defined by 

(M,0(r))(M',0'(r)) = (MM', 0(MV)0'(r)). 
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We denote the inverse image of F := SL2(Z) under the covering map by F := Mp2(Z). It 
is well known that F is generated by T := (( J } ) , 1), and S := (( ? ^ ) , V^)- 

Let iV be a positive integer. There is a certain representation p of F on C[Z/2NZ], the 
group ring of the finite cyclic group of order 2N. For a coset h G Z/2A^Z we denote by 
th the corresponding standard basis vector of C[Z/2NZ]. We write (-, ■) for the standard 
scalar product (antilinear in the second entry) such that {th, th') = Sh,h'- In terms of the 
generators T and S of F, the representation p is given by 

(2.1) p(T)(e,) 

(2.2) p{S){t,) 

Here the sum runs through the elements of Z/2A^Z and we have put e(a) = e^'^*". Note 
that p is the Weil representation associated to the one- dimensional positive definite lattice 
K = (Z, Nx^) in the sense of ^Bolj . |Brlj . |BruOj . It is unitary with respect to the standard 
scalar product. 

li k E |Z, we write Ml ^ for the space of C[Z/2A^Z]- valued weakly holomorphic modular 

forms of weight k for F with representation p. The subspaces of holomorphic modular 
forms and cusp forms are denoted by p and Sk^p, respectively. 

2.2. Harmonic weak Maass forms. In this subsection we assume that A; < 1. A twice 
continuously differentiable function / : EI — )■ C[Z/2A^Z] is called a harmonic weak Maass 
form (of weight k with respect to F and p) if it satisfies: 

(i) /(Mr) = 0(r)XM,0)/(r) for all (M, 0) G f; 

(ii) Afc/ = 0, 

(iii) there is a C[Z/2iVZ] -valued Fourier polynomial 

h {2N) neZ<o 

such that /(r) — P/(r) = 0(e~^^) as f — )■ oo for some e > 0. 
Here we have that 

is the usual weight k hyperbolic Laplace operator (see |BFlj ). The Fourier polynomial 
Pf is called the principal part of /. We denote the vector space of these harmonic weak 
Maass forms by Hk^p (it was called H^^ in [BFl] ). Any weakly holomorphic modular form 
is a harmonic weak Maass form. The Fourier expansion of any / G Hk^p gives a unique 
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decomposition f = + f , where 

(2.4a) /+(r)= c+(n, /^)gi^e,, 



h {2N) neZ 
n2>— oo 



(2.4b) /- {r)= ^~ (l - ^> 47r 



72 




4iV 





h (2N) raS 
n<0 



qiN e^. 



We refer to /"*" as the holomorphic part and to / as the non-holomorphic part of /. Note 
that c=^(n, h) = unless n = /i^ (4iV). 

Recall that there is an antilinear differential operator C, = C,k '■ H^ p — )• S2-k,p, defined by 



kdf 



(2.5) fiT)^af)ir):=2tv 

Here p denotes the complex conjugate of the representation p, which can be identified with 
the dual representation. The map ^ is surjective and its kernel is the space M\ p. There is 
a bilinear pairing between M2-k,p and Hk^p defined by the Petersson scalar product 

(2.6) {9j} = {9,m)--= [ {9,m)v'-'^'' 

Jr\m ^ 



for g G M2-k,p and / G If 9 has the Fourier expansion g = 'Yl,hn^i^i h)q^^^'^ th, and 

if we denote the Fourier expansion of / as in (12. 4p . then by |BFlt Proposition 3.5] we have 

(2.7) {g,f}= $^c+(n,/i)6(-n,/i). 

h {2N) n<0 

Hence {g, /} only depends on the principal part of /. 

2.3. The Shimura lift. Let k e lZ\Z. According to [EZl Chapter 5], the space Mk^p 
is isomorphic to Jk+i/2,N, the space of holomorphic Jacobi forms of weight A; + 1/2 and 
index N. According to [Ski] and |SZ] . Mk^p is isomorphic to J^'^^/2N^ ^^e space of skew 
holomorphic Jacobi forms of weight k + 1/2 and index A^. There is an extensive Hecke 
theory for Jacobi forms (see [EZj . |Skl] . [SZj). which gives rise to a Hecke theory on M^ p 
and Mfc p, and which is compatible with the Hecke theory on vector valued modular forms 
considered in |BrStj . In particular, there is an Atkin-Lehner theory for these spaces. 

The subspace S^^^ of newforms of Sk,p is isomorphic as a module over the Hecke algebra 
to the space of newforms S2^^\^{N) of weight 2k — 1 for ro(A^) on which the Fricke in- 
volution acts by multiplication with (—1)'^^^/^. The isomorphism is given by the Shimura 
correspondence. Similarly, the subspace SJ^^^ of newforms of Sk,p is isomorphic as a module 
over the Hecke algebra to the space of newforms S2^^\~{N) of weight 2k — 1 for ro(A^) on 
which the Fricke involution acts by multiplication with (— l)'=+i/2 (gee |SZj . [GKZj . |Skl] ). 
Observe that the Hecke L-series of any G G S2^^'-^{N) satisfies a functional equation under 
s ^ 2k — 1 — s with root number Eg = ±1- 

We now state the vector valued version of Theorem ll.li Let G G 5*2 '^"'(A^) be a nor- 
malized newform (in particular a common eigenform of all Hecke operators) of weight 2 
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and write Fq for the number field generated by the eigenvalues of G. If = — 1 we put 
p' = p, and if £:g = +1 we put p' = p. There is a newform g G SJ^J^p' mapping to G under 
the Shimura correspondence. It is well known that we may normalize g such that all its 
coefficients are contained in Fq. According to |BruO[ Lemma 7.3], there is a harmonic 
weak Maass form / G -f^i/2,p' whose principal part has coefficients in Fq with the property 
that 

Mf) = M-"9. 

This form is unique up to addition of a weakly holomorphic form in M[^^ ^, whose principal 
part has coefficients in Fq. 

In practice, the principal part of such an / can be computed as follows: We may complete 
the weight 3/2 form g to an orthogonal basis g , g2, ■ ■ ■ , gd of Sz/2,p' consisting of cusp forms 
with Fourier coefficients in Fq. Let / G -f^i/2,p' such that 

(2.8) {/, g} = 1, and {/, = for ^ = 2, . . . d. 

Then / has the required properties. In view of 02.71) the conditions of (12.81) translate into 
an inhomogeneous system of linear equations for the principal part of /. 

Theorem 2.1. Let G G 82^"^ (N) be a normalized newform. Let g G Sl^J^pn '^'^^ / ^ -^i/2,p' 
he as above. Denote the Fourier coefficients of f by c^{n,h) for n G Z and h G Z/2A^Z. 
Then the following are true: 

(1) If A ^ 1 is a fundamental discriminant and r G Z such that A = (mod 4A^) 
and > 0, then 

L(G, XA, 1) = 87r2||G|n|(?lpy® ■ c-(A)2. 

(2) If A 1 is a fundamental discriminant and r G Z such that A = (mod 4A^) 
and EqA < 0, then 

L'{G,XA,l) = ^ c+{-eGA,r)eQ ^ c+(-£gA, r) G F^. 

When 5'i/2,p' = {0} the above result also holds for A = 1, see also [Brum Remark 18]. 
This is for instance the case when is a prime. If is a prime and = —1, then the space 
Hi/2,p' can be identified with a space of scalar valued modular forms satisfying a Kohnen 
plus space condition. In that way one obtains Theorem 11.11 stated in the introduction. 

3. Computational aspects 

3.1. The automorphy method for vector valued weak Maass forms. To compute 
the Fourier coefficients of the harmonic weak Maass forms we use the so-called automorphy 
method, sometimes called "Hejhal's method". This is a general method which has been 
used to successfully compute various kinds of automorphic functions and forms on GL2 (M). 
It was originally developed by Hejhal in order to compute Maass cusp forms for the modular 
group and other Hecke triangle groups (cf. e.g. [He]). The method was later generalized by 
the second author in [Stlj to computations of Maass waveforms with non-trivial multiplier 
systems and arbitrary real weights, as well as to general subgroups of the modular group 
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(see also |St3j ). Another generalization to automorphic forms with singularities (Eisenstein 
series, Poincare series and Green's functions) was made by Avelin |AvH Rv2] . 

We will detail the adaptation of the algorithm to the case of vector-valued harmonic 
weak Maass forms for the Weil representation. 

For simphcity consider the representation p (the case of p is analogous) and A; G Z + |. 
Furthermore, in order to avoid questions of uniqueness we assume that either A; < or 
that k = ^ and that is prime. In these harmonic weak Maass form is uniquely 

determined by its principal part. For computational purposes it is not feasible to use the 
definition of p in terms of the action on the generators of the metaplectic group. We 
instead use formulas from [Stl] to evaluate p on the fixed (canonical) representative of 
M = (^^) G SL2(Z), i.e. p{M) := p{M,jM (r)) where jAf (r) = Vct + d is defined by the 
principal branch of the argument. 



3.1.1. The algorithm - phase 1. Let / G H^ p with a given (fixed) principal part Pf (r) = 
Pf,h (t) where Pf^h (t) = Y^n=-K ^ ^) (^^^ some finite A' > 0) and write 
f = f+ + f- (as in 2.3a and 2.3b) with /+ = Eh{27V) fh^h and f- = T.hi2N) fh^h where 







(r) = a (n, /i) g4jv -)- c"*" (n, /;,) g4jv and 

n=-K n>0 



n<0 



\AN 



for r = M + 2f G H. Our goal is to obtain numerical approximations to the coefficients 
c^(n,h). To formulate our algorithm we prefer to separate the u- and the f -dependence 
in / and therefore introduce the function W defined by W{v) = e"^'^'" if w > and 
W{v) = e"^™r(l - k,ATc\v\) if < 0. We also set c{n,h) = {n,h) for n > and 
c~ (ra, h) for n < and write e^^ (u) = . With this notation 







By standard inequalities for the incomplete gamma function one can show that 



1, v>0, 
(47r|t;|)"\ v<0, 



where Ck is an explicit constant only depending on k. To be able to determine a truncation 
point of the Fourier series above we also need bounds of the coefficients c{n,h). Using 
[BruFu, Lemma 3.4] it follows that there exists an explicit constant C > such that 

c{n,h) = O (exp (47rCv^)) , n +oo, 
c{n,h) = 0{\n\2), n — )■ — oo. 
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For A; < we are able to make the implied constants explicit using non-holomorphic 
Poincare series as in e.g. |Brlj or |He2j . For k = ^ we rely on numerical a posteriori tests 
to assure ourselves that the truncation point was choosen correctly. See e.g. Section 13.21 

Let e > and fix F < Iq = By the estimates above we can find an Mq = M {Y, e) 
such that the function / = J2h{2N) fh^h given by the truncated Fourier series 

A(r) = PM(r)+ c (n, (^) M 

0<|n|<Mo 



satisfies 



/>) - fir) 



< e 



for any T e Hy = {r en\Qr> Y}. Here \\z\\' = Y.h=i \zh\ for z E C^^. Let A 



c d , 



G SL2(Z) and set z = x + iy = At. Then y = ^At 



1 — —TTi < -T-^ < - and hence 



\3AiT)\ = \cT + d\ = I < p-. Using the fact that p is unitary it is now easy to see that if 
r, At G Hy then 



(3.1) 



/ (At) - J A {rf' p {A) f (r) < e (l + < 2e • 



2k 



Consider now a horocycle at height Y and a set of 2Q (with Q > Mq) equally spaced points 

1 - 2m 



4Q 



1-Q <m<Q. 



If we view the series fh as a finite Fourier series we can invert it over this horocycle and it 
is easy to see that if n is an integer with < \n\ < Mq and n = (4A^) then 



(3.2) ^ E A( 

^ m=l-Q 



W (^^) c {n, h) + a (n, h) e~^^ . 



One can also interpret the left-hand side as a Riemann-sum approximation to the integral 

fh (z) e^N {-nx) dx. 



Let = x*^ + iy*^ = T~^Zm (T^ G PSL2(Z)) denote the pull-back of z^ into the standard 
(closed) fundamental domain of PSL2(Z), = \^z = x + iy \ |x| < |, > l}. Using (13. ip 
we obtain 

fhizm)=jTjz*J J2 Pkh' (T^) fh' iz*J + l2eY-'% 

h' {2N) 

where phh' (^m) is the {h, h')-e\emeiat of the matrix p (T^), and we use |2ey~^'^] to denote 
a quantity bounded in absolute value by 2eY~'^''. Inserting this into (13. 2p we see that the 
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left-hand side can be written as 
Q 



2Q 



m=l- 



Phh' {Trr 

h' (27V) 



J2 h') exp (-j^y*mj em {lx*J 

l=-K ^ ^ 

+ c{l,h')w(^yC)e,r,{lx*J 
o<|;|<Mo ^ ^ 



eAN{-nXm) 



(3.3) 
where 



h' {2N) 0<|i|<Mo 



2feT 



= ^ (^m) (^m) W i^vCj emilxl, - nxm) and 

~ 1 ° f 2-111 \ 

^n-^H 5Z « X] (^-) (^-) ( "iiV^"^ ) ""'"^ ^^^"^ ~ ""^"^^ ■ 

^ fe' (2Ar) l=-K m=l~Q ^ ^ 

We thus have an inhomogcncous system of hncar equations which is (approximately) satis- 
fied by the coefficients c (n, h). Let T) = {{n^h) | < |n| < A'/q, < /i < 2A'"} (with a fixed 
ordering) and note that 12^1 = AMqN . If we set D = {d (n, /i))(„ fe)gxi' 

^ = ^ (^) = (Kf ) , , ^ , = - '^n^'^.^'W^ ({^y) and 

V / {h,n),{h',l)eV V4iV / 

W> = H> (F) = , < = W!: -a (n, /.) e'w^, 

we can write this linear system as I'D] linear equations in ll?! variables: 
(3.4) VD + W^6. 

In practice it turns out that the the matrix V is non-singular as soon as the subspace of 
p consisting of functions with a given singular part is one-dimensional. In these cases 
we can immediately obtain the solution as 

D = -V-^W, 

— * 

and since we know that the vector of the "true" coefficients, C — {c {n, h))^^ h)ev^ satisfies 



we see that 

C-D 



vc + w 



< \\v-'\ 



< 2eY 



-2k 



VC + W 



< 2eY 



-2k 



\v-'\ 



To obtain a theoretical error estimate we would thus need to estimate from below. 

Unfortunately this does not seem to be possible from the formulas above and we have to 
use numerical methods to estimate this norm. Hence, to obtain the Fourier coefficients up 
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to a (proven) desired precision we might have to go back and decrease the original e or 
increase either of Mq or Q. 

At this point one should also remark that the error bound is in general much 

worse than the actual apparent error, as verified by studying coefficients known to be 
integers. The reason for this is that the sums V^l"' exhibit massive cancellation and are 
therefore overpowered by the terms W (^^) on the diagonal. 



3.1.2. The algorithm - phase 2. Returning to (13. 3 p and solving for c (n, h) we see that 



(3.5) 



c (ra, h) 



-1 



h' {2N) \l\<Mo 



+ l2eT 



-2kl 



for any n, i.e. also when \n\ > Mq, provided that Q > M {Y). If we first choose Y such that 
W (^^) is not too small then we can in fact use this equation to compute c (n, h) with an 

error of size eW{^Y) . In this manner, we may produce long stretches of coefficients 
(before we need to decrease Y again) at arbitrary intervals Na <n< Nb without the need 
of computing intermediate coefficients above the initial set up to n = Mq. 

Remark 1. The exact same algorithm, with the non-holomorphic parts set to zero, also 
lets one compute holomorphic vector-valued modular forms for the Weil representation. 
This has been exploited by the second author, in verifying computations of holomorphic 
Poincare series in |RSS] . 

3.2. Heuristic error estimates. For A; < all implied constants and therefore all error 
estimates can be made explicit. In the remaining case which interests us. A; = |, the known 
bounds for the twisted Kloosterman sums are not enough to prove the necessary explicit 
bounds for the Fourier coefficients of the associated Poincare series. We are therefore not 
able to give effective theoretical error estimates in this case. However, this is not a serious 
problem since there are a number of tests we may perform on the resulting coefficients to 
assure ourselves of their accuracy. We list a few tests which we have used. 

• First of all, one can simply use two different values of Y and verify that the resulting 
vectors D = D{Y) are independent of Y. 

This test is completely general and can be used for all instances where the algorithm can be 
applied. Suppose now that we have a harmonic weak Maass form / G Hk^p of half-integral 
weight k such that (/) = with g G S2-k,p- We then know the following. 



• The coefficients -^/l A|c {—Eg ■ A) are proportional to the coefficients b{eG ■ ^) of g 
(cf. e.g. jBHIOl p. 3]). 

If additionally the Shimura lift of (7 is a newform G G S^^^j^ (ro(iV)) then we can predict 
that certain coefficients c"^(A) are algebraic (cf. e.g. |BruOt Sect. 7]) and if we are able to 
identify these coefficients as algebraic numbers to a certain precision this can be used as 
another measure of the accuracy. 
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3.3. Implementation. The first implementation of tlie above described algoritlim was 
made in Fortran 90, using tlie package ARPREC |ARj for arbitrary (fixed) precision aritfi- 
metic. Tlie second and more recent implementation was done in Sage |SA] . using the 
included package mpmath for arbitrary (fixed) precision arithmetic. The algorithms are 
currently under development but can be obtained on request from the authors. The final 
format we intend for these algorithms are standard classes for computing with vector and 
scalar-valued harmonic weak Maass forms in Sage or Purple Sage. 



4. Results 

4.1. Harmonic Maass forms corresponding to elliptic curves. In this section we 
present the numerical results we have obtained for harmonic weak Maass forms correspond- 
ing to weight two holomorphic forms associated to elliptic curves. We have concentrated 
on three particular examples. In Cremona's notation, these correspond to the curve llal 
of level 11 and the two curves 37al and 3761 of level 37. 

Recall that if the holomorphic weight 2 newform G of level has Atkin-Lehner eigen- 
value ±1 then the L-function L(G, s) has root number Eg = =Fl- Furthermore, since 
the root number of the twisted L-function L(G',xa,'S) is sign(A)xA(A^)^G and we always 
consider fundamental discriminants for which Xa(A^) = 1 we see that the central value 
L{G, xa, 1) vanishes if sign(A)£G = —1, i-e., if L{G, s) has an even functional equation we 
consider A < and otherwise A > 0. 

For each of these examples we computed a large set of central derivatives of the twisted 
L-functions with the appropriate A using Sage and the standard algorithms there which 
were developed by Dokchitser. We then fixed a harmonic weak Maass form with non-zero 
principal part Pf such that ,^3 (/) maps to G under the Shimura lift. In all cases we took 

a Poincare series P_a having principal part and computed an initial set of Fourier 

coefficients for this function using the methods described in the previous section. We then 
used the second phase of the algorithm and computed more Fourier coefficients. 

Note that for the results in this section, all initial "phase 1" computations were all 
performed using the new Sage package and all further, "phase 2" , computations were done 
in Fortran 90. 

We would like to give a flavour of the cpu-times involved. The initial computations, 
using our Sage code, took in all cases approximately 2 hours on a 2.66GHz Xeon processor. 
On the same processor, the cpu time for a single stretch of phase 2 calculations range 
between less than an hour for the smallest discriminant up to several days for the largest 
discriminant. 

As a measure of the accuracy of our computations one can consider the difference between 
the coefficients in Tables |21 Eland [Hand the nearest integer (the third column). To further 
support the correctness we also list, in Tables |31 IB) and normalized coefficients of the 
non-holomorphic parts, i.e. a/|A|c^(A)/ a/|Ao|c^(Ao) by some fixed non-zero coefficient 
of index Aq. 
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4.1.1. Hal. Here the unique newform of weight two and level 11 is given by 

G = 7]{Tf7]{llTf = q-2q^-q' + 2q^ + q' + ---e S^"" (ro(ll)) 

and the corresponding L-function L{G, s) has an even functional equation. Using Sage we 
computed all values of L'(G, xa, 1) for fundamental discriminants A < such that (^) = 1 
and |A| < 19703. This set consists of 2749 fundamental discriminants and amongst these 
we found 14 discriminants for which L'{G,xa,^) vanished up to the numerical precision 
(see Table E]). 

As a representative for the harmonic weak Maass form in the space -f^i/2,p corresponding 
to G, we choose the Poincare series P_5 with the principal part q~^{t7 — t-j). To compute 
the Fourier coefficients of P_5 we used the method described in the previous section with an 
initial e = 10~^° and Y = 0.5, which gave us a truncation point of Mq = 42, corresponding 
to A between —1847 and 1885. For a short selection of computed values of c"'"(A) see 
Table [Hand for a table of coefficients corresponding to all vanishing L'{G, xa, 1) see Table 
121 The first few normalized "negative" coefficients are displayed in Table El These values 
should be compared to the list in |Sk2| p. 505]. 

4.1.2. 37al. Consider the newform of weight two and level 37 which has an odd functional 
equation. The g-expansion is given by 

G = q-2q^-3q^ + 2q^ - 2q^ + - q^ + 6q^ + Aq^^ - 5g" + ■ • • G 5^"' (ro(37)) . 

Using Sage we computed all values of L'{G,xaA) ^^t fundamental discriminants A > 
such that (^) = 1 and |A| < 15000. This set consists of 2217 fundamental discriminants 
and amongst these we found 8 discriminants for which L'{G,xaA) vanished up to the 
numerical precision (see Table [5]). For the corresponding harmonic weak Maass form in 

3 

Hi/2,p we took P_3, which has a principal part q~'^{t2i + £21)- The initial computation 
was done in Sage, using e = 1 ■ 10~^^, which gave a value of Mq = 30, corresponding to 
discriminants in the range —4440 < A < 4585. For examples of the coefficients c"'"(A) see 
Tables H] and [5l The first few normalized "negative" coefficients are displayed in Table El 

4.1.3. 37hl. In this case we consider the newform of weight two and level 37 which has an 
even functional equation. The g-expansion is given by 

G = g + g3 _ _ g7 _ + 3^11 + . . . e ST"" (ro(37)) . 

Using Sage we computed all values of L'(G,xa,1) for fundamental discriminants A < 
such that (^) = 1 and |A| < 12000. This set consists of 1631 fundamental discriminants 
and amongst these we found 15 discriminants for which L'(G, xa,1) vanished up to the 
numerical precision (see Table [H]). For the corresponding harmonic weak Maass form in 
Hi/2,p we took P_i2, which has a principal part g~i48(e3o ~ ^so)- The initial computation 
was done in Sage, using e = 1 ■ lO^'^'^, which gave a value of Mq = 33, corresponding to 
discriminants in the range —4883 < A < 5029. For examples of the coefficients c"'"(A) see 
Tables [7] and [H The first few normalized "negative" coefficients are displayed in Table 1^741 
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4.2. Conclusions of the numerical experiments for weight two. In each of the ex- 
amples of weight two newforms that we studied we saw agreement with the theorem, i.e. 
the coefficients c"'"(A) (for fundamental discriminants with the appropriate property) were 
only algebraic when the corresponding central derivative L'{G,xaA) vanished. Further- 
more, we observed that in the cases we considered, the algebraic coefficients ^^(A) were in 
fact even rational integers. 

4.3. Further computations. To investigate whether a result analogous to Theorem 11.11 
also holds for newforms of weight 4, we computed L'{2, G, xa) for all newforms G of weight 

4 on ro(A^) whith 5 < < 150 and fundamental discriminants A with |A| < 300 and 
the property that the twisted L-function L{s, G, xa) has an odd functional equation. For 

5 < < 10 we additionally computed these values for fundamental discriminants A with 
|A| < 5000. Amongst all these values we did not find a single example of a vanishing 
derivative. Even though we did not get any positive case where we could test the theorem 
we still wanted to make sure that there was no easily accesible counter example. 

We therefore computed the Fourier coefficients, up to 40 digits precision, of the associated 
weight — I harmonic Maass form corresponding to all weight 4 newforms defined over Q 
for A^ up to 100. To test the accuracy (and making sure that the implementation was 
correct) we did not only rely on the provable error bounds, but also checked algebraicity 
of certain coefficients corresponding to non-fundamental discriminants. These coefficients 
were indeed all found to be integers or rational with fairly small denominators. In contrast 
to this, the Fourier coefficients corresponding to fundamental discriminants were found not 
to be similarly "simple" rational numbers. 

The L- value computations were performed in Sage |SA] . using the included version of 
Rubinstein's Icalc library [Ij. 

4.4. Tables. 

References 

[Avl] H. Avelin, Deformation of ro(5)-cusp forms, Math. Comp., 76 (2007), no. 257, 361-384. 
[Av2] H. Avelin, Computations of Automorphic Functions on Fuchsian Groups, PhD. Thesis, Uppsala 
University (2007). 

[AR] D. H. Bailey, Y. Hida, K. Jeyabalan, X. S. Li and B. Thompson, ARPREC: An arbitrary precision 

computation package. Version 2.2.4 (2009). 
[BSV] A. R. Booker, A. Strdmhergsson, and A. Venkatesh, Effective computation of Maass cusp forms. 

Int. Math. Res. Not. ID71281 (2006), 1-34. 
[Bol] R. Borcherds, Automorphic forms with singularities on Grassmannians, Invent. Math. 132 (1998), 

491-562. 

[BOl] K. Bringmann and K. Ono, The f{q) mock theta function conjecture and partition ranks. Invent. 

Math. 165 (2006), 243-266. 
[B02] K. Bringmann and K. Ono, Arithmetic properties of coefficients of half-integral weight Maass- 

Poincare series. Math. Ann. 337 (2007), 591-612. 
[B03] K. Bringmann and K. Ono, Dyson's ranks and Maass forms. Annals of Math. 171 (2010), 419-449. 
[Brl] J. H. Bruinier, Borcherds products on 0(2, 1) and Chern classes of Heegner divisors, Springer 

Lecture Notes in Mathematics 1780, Springer- Verlag (2002). 
[Br2] J. H. Bruinier, Harmonic Maass forms and periods, in preparation. 



Table 1. E = Hal, P_5 e Hi 



-(A) 



-7 2.8463370190285980186651576519711751393948073861988 

-8 2.5138482002575729165892711124435774460030012341762 

-19 4.8192428148963255861870924437043119000673519519042 

-24 5.1018494088339703187507177747597327598232414476402 

-35 4.7515892636101723769649079639675162004017245362399 

-39 -1.6466690697010481166272091028219327677356442914804 

-40 1.1470941388138074683747768314723689292860539962900 

-43 -1.7622439638503327722737780360046423237568367805048 

-51 2.0736222999878741718629718432682995552582880786065 

-52 1.5723528683914990387103216700146317562411438497615 

-68 9.6889673322938493992006043404127469979247370067926 

-79 1.7557351755436160739388564340027760291317089229254 

-83 -7.1767664383427675609861242907417950544611683162859 

-84 6.1666200626587315159968126799603650525586539365601 

-87 -7.7230036424433334541697484050338439023979647483280 

-95 7.8467572084064151556661839046504144426199227897994 

-811 3.0046247983067285336553431175489765847382042907105 

-820 -6.0493754250387304262091147332158578046510749019315 

-824 -5.7985199999999999999999999999999999999999999999999 

-827 1.8535489407871222859528059423736067736521222528554 

-831 -6.7911392225835416131083026699411310608420151994986 

-996 -3.5516294505685820400211045047063422129168082892941 

-1003 1.0811934742079303073802766406181476437608668928409 

-1007 -3.9469248000000000000000000000000000000000000000000 

-1011 -3.7685140824429636488934775010060106547341275101054 

-1019 3.3790315957549749442218769650593817997888628818338 



i^'(G,XA,l) 


.1000 


1.22556687406888 


1000 


• 10°° 


1.88791354720393 


looo 


■ 10°° 


7.51391655667451 


lOOO 


■ 10°° 


4.02744559024300 


lOOO 


.1000 


7.64786334637073 


1000 


■1001 


2.97721567216550 


lOOO 


■10°i 


5.58789208952436 


lOOO 


•lOOl 


1.18814465355690 


lOOi 


•1001 


1.30416363302768 


1001 


•1001 


5.14853759817659 


1000 


.1000 


3.80344864881298 


lOOO 


•lOOi 


4.75620653690677 


lOOO 


.1001 


6.43843846621214 


1000 


•1001 


6.53746327159376 


lOOo 


■ lOOi 


2.35584785481347 


lOOO 


■ 1001 


3.03660486030085 


lOOO 


• 1006 


1.25949136911120 


lOOi 


.1006 


1.19119437485937 


1001 


■ 1006 


-6.6 ■ 10-24 




■ 1006 


1.60961273159300 


lOOi 


•1006 


3.36744068632019 


lOOO 


•1007 


1.15828152096335 


1001 


• 1007 


2.53076681967579 


1001 


■1007 


-1.1 ■ 10-22 




•1007 


1.84592490209627 


lOOi 


■ 1007 


1.68145450009782 


lOOi 
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Table 2. E ^ Hal, P_5 e Hi ^ 



A 


C {IS) 


\r+(A) 
\C [lA) 


— \r 


V 


—824 


—5798520 


3.0 


• 10" 


-76 


—1799 


—2708450784 


2.7 


• 10 


-46 


—4399 


—68135748249936640 


2.3 


■ 10 


-21 


-8483 


214445760716391388216704 


9.1 


■ 10" 


-28 


-11567 


-12412267149099919205092899456 


1.6 


■ 10" 


-ZO 


-14791 


66850179291021019012709832099520 


3.1 


• 10- 


-30 


-15487 


-478732239405182448762415030881280 


5.6 


• 10- 


-30 


-15659 


-804489814454597618648064770159415 


6.8 


■ 10- 


-30 


-15839 


-1162122495004344641799524116135680 


7.0 


■ 10- 


-30 


-16463 


4542575922533728228643934862230144 


1.5 


• 10- 


-30 


-17023 


-23302350713109514450879400185948800 


2.0 


• 10- 


-29 


-17927 


110133238181959291703634158808374784 


1.2 


• 10- 


-29 


-18543 


464726791864282489334104058164482624 


1.9 


• 10- 


-29 



Table 3. E = Hal, P_5 e .ffi^^. Coefficients are scaled by c (1). 



A 


VAc-(A) 


|c-(A) 


-[c-(A)]| 


4 


-3 


2.0 


IQ-lOO 


5 


5 


2.1 


10-99 


9 


-2 


1.7 


10-100 


12 


5 


8.0 


10-100 


16 


4 


1.5 


10-99 


20 


5 


1.1 


10-100 


25 





1.0 


10-100 


36 


6 


1.0 


10-99 


37 


5 


4.2 


10-99 


45 





6.4 


10-99 
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Table 4. E ^ 37al, P_3 e Hi „ 

2 'P 



(A) 



L'{G,XA,1) 



1 -2.8176178498959956879756075537515493438922975370716 

12 -4.8852723826201225228227029607337071669095284814788 

21 -1.7273925723265275652082007397068992218426924398791 

28 6.7819399530394779828450578400669420938246859928076 

33 5.6630232015906998168220545669245622604190884430064 

37 -9.1326561374611652958506448407204050631184401026129 

40 4.0098509269543637915254766073122850557290259963615 

41 6.5637495744757231959699415722023547525400239084778 
44 9.6886404434506397321859573425794267139455920322171 
53 -5.6688852568232517859984506645723944339238503996096 
65 -6.0328072889521477971996798071175156059595671497733 
73 3.4874711835362408853804154923777452565842552803223 
77 2.2699132373705254600799448087564660809534768699467 
85 -7.6894617048676272061865441758881289699552927122551 

1481 -3.2715595098273932057423414526408419506801164996884 

1484 -1.3432792297590353562651264178555980321660674399890 

1489 8.9999999999999999999999999999999999999999999999999 

1496 1.1199440423162819213593329208218112792285448658029 

1501 -5.8188238119388864901078937905792951783427273426771 

4376 -3.6731327299348159991042234350611468700145535059868 

4377 -5.0062522276143084997015960658866819832113068397294 
4393 6.6000000000000000000000000000000000000000000001468 
4396 -2.3023069110811173762943326075771836710063196221488 
4412 -3.1500483730098996665306117169085504925545562420809 



10-01 

10-01 

10-01 
10-01 

10-01 

10-01 

10-01 
10-01 

10-01 

10-01 

10-01 

10-01 
10-01 

10-01 
10"o 

lOOi 
1000 

lOOi 

1002 
1002 
1002 

1001 

1002 
1002 



3.05999773834052 

4.29861479867736 

9.00238680032537 

4.32726024966011 

3.62195679113882 

3.47328771649229 

3.70588717878444 

5.93680171871573 

1.01334656625280 

2.61746665637296 

7.67818286326206 

2.92507284795068 

3.42067600398534 

9.90133670369251 

5.26994449124823 

3.86746474997364 

-3.7 • 10-23 

2.27616829409607 

6.06007663972706 

2.03155740209437 

2.27150950159608 

5.8 • 10-23 

2.00437958330233 

3.73011222569745 



10-01 

lOOO 

lOOO 

1000 

1000 

lOOO 

lOOO 

1000 

1001 

lOOi 
lOOO 

lOOO 

1010 

lOOO 

lOOO 
lOOi 

lOOi 
lOOO 

lOOi 
lOOO 

lOOO 
lOOi 
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Table 5. 37al, P_3 e Hi „ 

2 'f^ 



A 
ZA 


c (Aj 


\c (Aj 






1489 


9 


1.6 


10" 


-72 


4393 


66 


1.5 


10- 


-45 


5116 


-746 


8.5 


10- 


-23 


5281 


153 


8.2 


10" 


-23 


5560 


-1124 


1.2 


10- 


-22 


5761 


-974 


1.1 


10- 


-22 


6040 


-1404 


4.2 


10- 


-23 


6169 


336 


1.1 


10- 


-22 



Table 6. E ^ 37al, P_3 e Hi p. Coefficients are scaled by VSc~{-3). 

A V\A\c-i^) |c-(A)-[c-(A)]| 



-4 


1 


4.0 


lQ-84 


-7 


-1 


5.0 


10-84 


-11 


1 


4.5 


10-84 


-12 


-1 


2.0 


10-84 


-16 


-2 


1.1 


10-83 


-27 


-3 


1.3 


10-83 


-28 


3 


1.4 


10-83 


-36 


-2 


1.0 


10-83 


-40 


2 


3.6 


10-85 


-44 


-1 


1.1 


10-83 
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Table 7. E = 3761, P_i2 e Hi 



(A) 



-3 1.0267149116920354474451980654263083626994977508118 

-4 1.2205364009670316625279102409757685190938519711297 

-7 1.6900297463200076214148752932012403965170838158011 
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Table 8. E ^ 37W, P_i2 e Hi ^ 
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Table 9. = 37&1, P_i2 e Hi p. Coefficients are scaled by c"(l). 
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